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ABSTRACT 

We present measurements of anisotropy in the Cosmic Microwave Background (CMB) from the first season 
of observations with the Degree Angular Scale Interferometer (DASI). The instrument was deployed at the South 
Pole in the austral summer 1999-2000, and made observations throughout the following austral winter We present 
a measurement of the CMB angular power spectrum in the range 100 < Z < 900 in nine bands with fractional 
uncertainties in the range 10-20% and dominated by sample variance. In this paper we review the formalism used 
in the analysis, in particular the use of constraint matrices to project out contaminants such as ground and point 
source signals and to test for correlations with diffuse foreground templates. We find no evidence of foregrounds 
other than point sources in the data, and find a maximum likelihood temperature spectral index (3 = -0.1 ±0.2 
(1 cr), consistent with CMB. We detect a first peak in the power spectrum at / ^ 200, in agreement with previous 
experiments. In addition, we detect a peak in the power spectrum at / ^ 550 and power of similar magnitude 
at I ^ 800 which are consistent with the second and third harmonic peaks predicted by adiabatic inflationary 
cosmological models. 
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1. INTRODUCTION 

Subtle temperature fluctuations in the Cosmic Microwave 
Background (CMB) radiation, first observed by the COBE 
DMR experiment (Smoot et al. 1992), offer a glimpse of the 
early Universe at the epoch of matter-radiation decoupling, well 
before the non-linear gravitational collapse of matter led to the 
structure we see in the Universe today. The theory of inflation, 
originally proposed as an explanation for apparent isotropy in 
the CMB and a spatial curvature near unity (Guth 1981; Linde 
1982; Albrecht & Steinhardt 1982), was further developed to 
predict nearly scale-invariant adiabatic density perturbations 
arising from quantum fluctuations as the source of inhomogene- 
ity in the Universe (Guth & Pi 1982; Hawking 1982; Starobin- 
sky 1982; Bardeen, Steinhardt, & Turner 1983). On angular 
scales causally disconnected at the epoch of decoupling (> 1°), 
anisotropy in the CMB reflects the primordial inhomogeneity of 
matter (see, e.g.. White, Scott, & Silk 1994). On scales smaller 
than the sound horizon at the time of decoupling, the primordial 
adiabatic density perturbations force gravity-driven acoustic os- 
cillations in the photon-baryon fluid which lead to a harmonic 
series of peaks in the CMB angular power spectrum (Bond & 
Efstathiou 1984; Vittorio & Silk 1984). A detection of a series 
of peaks in the CMB angular power spectrum would provide 
strong evidence for the inflationary view of the early Universe, 
and would make alternate theories of structure formation such 
as defect models difficult to support (Albrecht et al. 1996). A 
third harmonic peak of comparable or greater magnitude than 
the second is the signature of baryonic matter re-compressing in 
dark matter potential wells — a detection of significant power 



in the third peak region would strongly support the presence of 
dark matter in the Universe (Hu & White 1996). Within the 
context of standard cosmological models, the observed CMB 
angular power spectrum can be used to determine fundamental 
cosmological parameters (Knox 1995; Jungman et al. 1996). 

The angular scale of the first acoustic peak, coupled with 
the sound horizon size at last scattering, provides a classical 
angular-diameter distance measure of spatial curvature. Recent 
experiments have determined the location of the degree-scale 
first peak, providing strong evidence for a spatially flat Universe 
(Mifler et al. 1999; Mauskopf et al. 2000; de Bernai-dis et al. 
2000; Hanany et al. 2000). Additional parameters such as the 
baryonic matter density (ilt, h^) can be extracted by resolving 
the higher-order peaks in the power spectrum. Hints of struc- 
ture appear to be present in some smaller angular scale CMB 
power spectra (Hanany et al. 2000; Padin et al. 2001a), but 
there has been no clear detection of the predicted higher-order 
acoustic peaks. 

Until recently, most observations of degree-scale CMB 
anisotropy have been made with ground and balloon based 
single-dish experiments. The Degree Angular Scale Interfer- 
ometer (DASI), like its sister instrument the CBI (Pearson et al. 
2000) and the VSA (Jones 1996), is a new compact interfer- 
ometer constructed specifically for observations of the CMB. 
Since interferometers are inherently differencing instruments, 
they allow measurements of anisotropy in the CMB that are 
free of many of the systematics which need to be controlled 
in non-interferometric CMB experiments. DASI is designed to 
measure the angular power spectrum of the CMB and produce 
high signal-to-noise images on angular scales corresponding to 
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the first three acoustic peaks predicted by adiabatic inflationary 
models for a flat Universe. 

In this, the second of three papers describing the results of 
the first season of DASI observations, we focus on the deter- 
mination of the CMB power spectrum from the calibrated data. 
This paper includes a discussion of the data analysis method, 
potential sources of astronomical and terrestrial contamination 
in the data, tests of data consistency, and the resulting angular 
power spectrum. Detailed descriptions of the instrument, ob- 
servations, and data calibration are given in Leitch et al. (2001, 
hereafter Paper I). The extraction of the cosmological param- 
eters from the DASI angular power spectrum is presented in a 
third paper, Pryke et al. (2001, hereafter Paper III). 

2. INSTRUMENT 

An extensive discussion of the instrument design and its per- 
formance are given in Paper I; here we provide a brief overview 
stressing the aspects of the instrument that are particularly rele- 
vant to the measurement of the CMB angular power spectrum. 

Interferometers offer several unique features which make 
them well suited for sensitive measurements of CMB 
anisotropy: 1) the Fourier transform of the sky plane is mea- 
sured directly; the effective sky brightness differencing is in- 
stantaneous, 2) 180° phase switching at the receivers with syn- 
chronized post-detection demodulation is used to reduce instru- 
mental offsets (DASI uses both fast (40 KHz) phase switching 
with hardware demodulation and slow (^ 1 Hz) phase switch- 
ing with software demodulation to reduce offsets to well below 
the /xK level), and 3) the effective window function in /-space is 
well understood and uncertainties in the primary beam response 
do not lead to uncertainties in the resultant power spectrum that 
become large at small angular scales. 

An interferometer directly samples the Fourier transform of 
the sky brightness distribution. The response pattern on the sky 
for a given pair of antennas is a sinusoidal fringe pattern at- 
tenuated by the primary beam of the individual antennas. For 
a pair of antennas with a physical separation (baseline) vector 
b, the center of the measured Fourier wavevector components, 
labeled u or (m,v), is given by m = fo^/A and v = by/X, where 
bx and b^. are the projections of the baseline normal to the line 
of sight, and A is the observing wavelength. The approximate 
conversion to multipole moment is given by / « 27r |u| (White 
et al. 1999a), with a width A/ that is related to the diameter of 
the apertures in units of the observing wavelength (see §5.1). 

Both real and imaginary components of the Fourier plane can 
be measured with a complex correlator. The real (even) compo- 
nent is measured by correlating the pair of signals without a rel- 
ative phase delay; the imaginary (odd) component is measured 
by correlating the signals with a 90° phase shift introduced into 
one of the signal paths. The averaged correlated output of the 
interferometer is called the visibility (see eq. [1]) and is the fun- 
damental data product. 

DASI is a 13-element interferometer operating at 26- 
36 GHz. The 13 antenna elements are arranged in a threefold 
symmetric configuration on a common mount which can point 
in azimuth and elevation. The mount is also able to rotate the 
array of horns about the line of sight to provide additional (m, v) 
coverage as well as the ability to perform consistency checks. 
The 13 elements provide 78 baselines with baseline lengths in 
the range 25-121 cm. The configuration of the horns was cho- 
sen to provide dense coverage of the CMB angular power spec- 
trum from 100 </< 900. 



Each DASI antenna consists of a 20-cm aperture-diameter 
lensed corrugated horn which defines the ^ 3.°4 FWHM field of 
view of the instrument. The receivers use cooled low-noise high 
electron mobility transistor (HEMT) amplifiers (Pospiezalski 
et al. 1995), and have system noise temperatures referred to 
above the atmosphere ranging from 18 K to 35 K at the center 
of the band. The receivers downconvert the 26-36 GHz RF 
band to a 2-12 GHz IF band. Each receiver IF is further spUt 
and downconverted to ten 1 GHz wide bands centered at 1.5 
GHz. An analog correlator (Padin et al. 2001b) processes the 1 
GHz bands into 780 complex visibihties. 

The stability of the instrument, its location at the South Pole, 
and the fact that its mount is fully steerable, have given us great 
flexibility in designing and adapting our observing strategy. We 
are able to choose fields to avoid foreground contamination, 
balance sensitivity and sky coverage, and observe in patterns 
that reject ground and other spurious signals while produc- 
ing datasets containing correlations which are computationally 
tractable. 

3. OBSERVATIONS 

CMB fields were observed during the period spanning 05 
May-07 November 2000. The data presented here comprise 
97 days of observation, representing an observing efficiency of 
better than 85% (of the days devoted exclusively to CMB ob- 
servations), with the remainder lost to hardware maintenance 
and repairs. Observations were never prevented due to weather, 
and only 5% of data were lost due to weather based edits, con- 
firming previous assessments of the exceptional quality of the 
site (Lay & Halverson 2000; Chamberlin, Lane, & Stark 1997). 

The presence of near-field ground contamination at a level 
well above the CMB signal limits our ability to track single 
fields over a wide range in azimuth. Repeated tracks over the 
full azimuth range show a strong variation of the ground with 
direction, with amplitudes of tens of Jy on some of the shortest 
baselines, but there is little evidence for time variability on pe- 
riods as long as five days. One advantage of observing near the 
South Pole is that sources track at a constant elevation, which 
enables us to observe several sources at constant elevation over 
a given range in azimuth. Observations were divided among 4 
constant declination (elevation) rows of 8 fields, on a regular 
hexagonal grid spaced by Ih in right ascension, and 6° in decU- 
nation. The grid center was selected to avoid the Galactic plane 
and to coincide with a global minimum in the IRAS 100 /mum 
map of the southern sky. Each field in a row was observed over 
the same azimuth range, leading to a nearly identical ground 
contribution. The elevation of the rows are 61°, 67°, 55°, 49°, 
which we label the A, B, C and D rows for the order in which 
they were observed (see Paper I for full coordinates). The field 
separation of Ih in RA represents a compromise between im- 
munity to time variability of the ground signal and a desire to 
minimize inter-field correlations. 

A given field row was observed daily over two azimuth 
ranges, for a total of 16 hours per day, with the remainder of the 
time divided among various calibration and pointing tasks (see 
Paper I). Phase and amplitude calibration were accomplished 
through observations of bright Galactic sources, permitting de- 
termination of the calibrator flux on all baselines to better than 
2%. Absolute pointing error determined by offsets between 
DASI detected point source positions and PMN southern cat- 
alog coordinates (Wright et al. 1994) was less than 2', with a 
drift <C 1' over the period during which each row was observed. 



3 



The number of days for which each of the four rows was ob- 
served is 14, 24, 28 and 31 for the A, B, C and D rows, respec- 
tively, for a total integration time of 28-62 hours per field. 

4. CALIBRATION & DATA REDUCTION 

Absolute calibration of the telescope was achieved through 
measurements of external thermal loads; the cahbrations were 
then transferred to bright astronomical sources. The flux 
scales resulting from two independent calibrations performed in 
February 2000 and February 2001 are found to agree to 0.3%, 
consistent with our estimate of 1% overall statistical uncertainty 
in the measurement and transfer procedure. The systematic un- 
certainty in determining load coupUng and effective tempera- 
ture is 3%, which is the dominant contribution to the uncer- 
tainty in our overall flux scale. This uncertainty, expressed as a 
percentage of Q, is 7% at 1 cr and is constant across all power 
spectrum bands. Band-power measurements are also affected, 
though weakly, by errors in the estimated aperture efficiency, 
on which our uncertainty is 4% (see Paper 1). This contributes 
a band-power uncertainty which is constant at 4% except in the 
three lowest-/ bands, where a cancellation of errors causes it 
to decrease. In using the current DASI results for parameter 
estimation (Paper 111), we have found no significant difference 
between treating this small beam uncertainty separately with 
its low-/ variation included, and folding it together with the /- 
independent flux scale uncertainty. We therefore adopt a total 
combined calibration uncertainty of 8% (1 cr), expressed as a 
fractional uncertainty on the C/ band powers (4% in AT /T). 

Raw data from the correlators, along with monitoring data 
from various telescope systems, are accumulated in 8.4-s inte- 
grations. These short integrations are edited before being com- 
bined for analysis. Baselines are rejected for which the phase 
offset or relative gain between the real and imaginary multi- 
pliers exceed nominal values. Data are also rejected when an 
LO has lost phase lock, when a receiver has warmed, or to trim 
field scans so that all eight fields are observed over precisely 
the same azimuth range. We also edit data for which noise cor- 
relations between baselines indicate strong atmospheric fluctu- 
ations. 

The edited and caUbrated data are combined into 1-hr bins, 
with uncertainty in the bins estimated from the sample variance 
of the 8.4-s integrations. In order to implement ground con- 
tamination common mode rejection, it is necessary that a given 
visibility be measured for all 8 flelds in a row; we cut all base- 
lines that do not satisfy this criterion. We apply more stringent 
edits for (m,v) radii < 40, which we find are more susceptible 
to contamination. For these visibilities, we retain only data for 
which both the sun and moon are below the horizon. To mini- 
mize the risk of biasing the power spectrum results, we do not 
edit the data based on the level of the signal. We have varied the 
threshold values of the weather, calibrator, and lunar/solar edit 
criteria with no significant effect on the results. Collectively, 
these edits reject about 40% of the data. See Paper I for a more 
comprehensive description of the data edits. 

All observations of a given set of fields are then combined, 
and it is these 1560 combined visibilities per field (78 complex 
baselines x 10 correlator channels, before edits) which form 
the input to the angular power spectrum UkeUhood analysis. 

5. ANALYSIS 
5.1. Formalism 



The DASI instrument makes direct measurements of the 
Fourier plane, and the angular power spectrum can be extracted 
from the data without creating an image. The caUbrated output 
of the interferometer is the visibility. 



2kBT ~ AT 
V(u)=^g(X)A(u,X)*—(u), 



(1) 



which is the convolution of the Fourier Transform of the sky 
brightness distribution, AT(u)/T, with the antenna aperture 
field autocorrelation function, A(u, A), and ^(A) is a 2% cor- 
rection between the Rayleigh- Jeans and Planck functions. The 
aperture field autocorrelation function A(u, A) is radially sym- 
metric, peaking at [u| = 0, and tapering smoothly to zero at 
|u| = D/X, where D is the aperture diameter. In the flat-sky 
limit, which is appropriate for the 3.°4 FWHM DASI flelds, 

2v 

AT \ 

' ' ~ q for lul > 10 



5(|u|) = 



i=27r|u| 



(2) 



(White et al. 1999a; Hobson, Lasenby, & Jones 1995); we as- 
sume / = 27r |u| over the /-range to which DASI is sensitive. For 
a single visibility, a simple quadratic estimator <S of the quantity 
2w |u|^5(|u|) » lil + l)Ci/i2TT) is given by (White et al. 1999b) 
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<S = 27r|u|- 
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(4) 



where is the instrument noise variance for the measured vis- 
ibility y(u) and equation (4) gives the number speciflc to the 
DASI apertures and a power spectrum in units of /zK^. The vari- 
ance of the visibility is thus directly related to C/ centered at the 
basehne length |u| = l/(2n), with width A|u| ~ 12 (FWHM) 
determined by the width of the aperture fleld autocorrelation 
function. 

While the simple quadratic estimator above is useful for un- 
derstanding the relationship between the visibility and the an- 
gular power spectrum, we have chosen a maximum hkehhood 
method in the present analysis. We have adopted the iterated 
quadratic estimator approach of Bond, Jaffe, & Knox (1998) 
to find the maximum likelihood values of the angular power 
spectrum for a piecewise flat /(/+ l)C;/(27r) power spectrum in 
nine bands. A data vector A of length A^ = 1560x32 (before 
data edits) is constructed by combining observations of each 
visibility for each of the 32 flelds. The UkeUhood function for 
a set of parameters k is 

^^^^^ = ;;;;T a72 exp(-i A^C(K)-i A), (5) 

(27r)'v/2|C(«;)|''^ 

where the covariance matrix 

C{k) = Ct{k)+C„+Cc (6) 
is the sum of the theory, noise, and constraint covariance ma- 
trices, described below, and is a function of the parameters. 
The parameters k which we estimate are the band powers, 
2iTU^S{u) w /(/ + l)Q/(27r). The theory covariance matrix, Cj, 
is given by 

Cr(At) = (VV^) = ^Atp5p, (7) 

p 

where V is the vector of noiseless theoretical visibilities and the 
sum is over the piecewise flat bands p. The matrices Bp repre- 
sent the instrument fllter functions to fluctuation power on the 
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sky. They are constructed from the overlap integral of the aper- 
ture field autocorrelation functions of pairs of baselines where 
they sample the same Fourier modes on the sky, 

1 {2kBTf 



27r 



1 /■!« 



I»2|, ^ 

L |uf 



A(u,-u',A,)[A(u,-u',A,) 



±A(u,+u',A,)]. (8) 
Here |ui |^ , \u2\p are the lower and upper radial limits of band 
p, and + and - are used for the real and imaginary parts of 
the visibiUty, respectively, which we treat as separate elements 
in the data vector. The vectors u, and Uj are the (u,v) coor- 
dinates of the visibility data vector elements A,- and Aj re- 
spectively (White et al. 1999a, 1999b; Hobson et al. 1995). 
Equation (8) is simply an overlap integral between an aper- 
ture autocorrelation function centered at u, and ones centered 
at iuy; only one of the two overlap terms is non-zero, except 
for the shortest baselines, and the total integral is non-zero only 
when u, ± < D/Xi+D/Xj, i.e., when the two visibilities are 
sensitive to some of the same wavevector components in the 
(m,v) plane. We use the theoretical aperture fields in this cal- 
culation, which is justified by the good agreement between the 
theoretical and measured beams (see Paper I and Halverson & 
Carlstrom 2001). Fields are separated such that the inter-field 
data vector elements are essentially uncorrelated except in the 
highest elevation row, for which we calculate the appropriate 
correlations. 

The noise covariance matrix of the combined visibihty data 
vector is diagonal, with elements C„„ = afj estimated from the 
sample variance in the 8.4-s integrations over the 1-hr obser- 
vations. To verify the assumption that C„ is diagonal, we have 
calculated the sample covariance matrix from the data in each 
of the ten frequency channels for all 1-hr observations. We find 
rare occasions where the visibiUties are strongly correlated due 
to atmospheric fluctuations. Our weather edits consist of cuts 
based on the strength of these correlations; we cut observa- 
tions in which any off-diagonal correlation coefficient exceeds 
±0.36, but the data consistency does not depend strongly on 
this value. 

To reduce near-field ground contamination and point source 
contributions to the power spectrum, we employ the constraint 
matrix formahsm described in Bond et al. (1998) to marginalize 
over potentially contanoinated modes in the data. Specifically, 
for a given mode q, we construct a constraint matrix from the 
outer product of a template vector Tq, 

Cc« = TqT/ (9) 

and 

Cc = aYCca, (10) 



where a is a number large enough to de-weight the undesired 
modes without causing the covariance matrix to become singu- 
lar to working precision. In practice we have found that a can 
range over at least seven orders of magnitude without affecting 
the results. As an example of a template vector, in the sub- space 
of the data vector consisting of a single visibility observed in 
eight fields, a template vector Tq = [11111111]^ is used to con- 
strain a common mode with the same amphtude in all eight 
fields during any given 8-hr period of observation. This effec- 
tively rejects a constant amplitude component, even one which 



varies in amplitude between subsequent 8-hr observations, such 
as a temporally drifting noise component with a period of many 
days. Any mode in the data which can be described as a rela- 
tive amplitude between data vector elements, as in the exam- 
ple above, can be constrained. We use this method to reduce 
near-field ground contamination in the field rows (see §5.2) 
and contributions from point sources with known positions (see 
§5.3). It is equivalent to marginalizing over these modes, with 
no knowledge of their amplitude scale. For each point source 
we also use constraint matrices to marginalize over arbitrary 
spectral indices, which we approximate as an amphtude slope 
across the ten frequency channels. We emphasize that we do 
not subtract ground components or point sources from the data. 
Instead we render the analysis insensitive to these modes in the 
data using the methods described above. 

The covariance matrix is block diagonal by field row, which 
allows us to invert the four sub-matrices in parallel. We further 
compress the matrix by combining visibilities and covariance 
matrix elements from adjacent frequency channels, which are 
nearby in the (m,v) plane and are therefore highly correlated. 
This data compression has the effect of slightly increasing the 
uncertainties and the anti-correlations between the resulting 
band-power estimates, but preserves the original piecewise-flat 
theoretical power spectrum model. 

The UkeUhood analysis software was extensively tested 
through analysis of simulated data. The analysis software and 
data simulation software were written by different authors in 
order to check for potential errors in the analysis. Omitting the 
constraint matrix leaves a sparse covariance matrix which can 
be rapidly inverted, and we can analyze a simulated data set in 
a few minutes of CPU time. Multiple simulated datasets were 
generated from an input model power spectrum, each with in- 
dependent sky and instrument noise realizations; the analysis 
software recovered the input model power spectrum within the 
estimated uncertainties, and these uncertainties were found to 
match, on average, the scatter in the band-power estimates over 
many realizations of the simulated data vector. We found no 
evidence of bias in the maximum likelihood band-power esti- 
mators. Ground signal and point source constraints were tested 
by constraining these modes in simulated data which contained 
both ground and point source components; both components 
are effectively eliminated, and the constraint matrices do not 
introduce artifacts into the power spectrum. 

5.2. Ground Constraints 

To remove sensitivity to the ground signal, we apply a con- 
straint which marginalizes over a common component across 
eight fields for each visibility, as described above. Addition- 
ally, using sensitive consistency tests described in §6.1, we find 
evidence of a temporally drifting component of the ground sig- 
nal on 1- to 8-hr time scales, subtle but present for all baselines 
and noticeably stronger for short baselines. We therefore ap- 
ply a linear drift constraint to all visibilities, and a quadratic 
constraint for |u| < 40. The additional constraints have little 
effect on the power spectrum, which makes us confident that 
sensitivity to ground signal is effectively eliminated. 

5.3. Point Source Constraints 

As predicted for our experimental configuration (Tegmark & 
Efstathiou 1996), point sources are the dominant foreground 
in the DASI data. To remove point source flux contributions 
using the constraint matrix formalism above, we require only 
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the positions of the sources, not their flux densities. We con- 
strain 28 point sources detected in the DASI data itself, in which 
we can detect an 40 mJy source at beam center with > 4.5 cr 
significance. The estimated point source flux densities range 
from 80 mJy to 7.0 Jy. We also constrain point sources from 
the PMN southern (PMNS) catalog (Wright et al. 1994) with 
4.85 GHz flux densities, ^5, which exceed 50 mJy when multi- 
plied by the DASI primary beam. We use this flux density Umit 
for the constrained point sources because the loss of degrees of 
freedom resulting from the inclusion of all point sources in the 
PMNS catalog would be prohibitively large. We have tested for 
the effect of possible absolute pointing error by displacing the 
point source position templates. A uniform displacement of the 
PMNS catalog coordinates by less than our estimated pointing 
error of 2' (see Paper I) does not have a significant effect on 
the angular power spectrum, except in the three highest-/ bins 
where the effect is ~ 10%. For the brightest point sources, po- 
sitions accurate to < 1' are required. We can extract positions 
to the necessary accuracy from the DASI data (see Paper I). 

In addition to the point sources constrained above, we make 
a statistical correction for residual point sources which are too 
faint to be detected by DASI or included in our PMN source ta- 
ble. To do this, we estimate the point source number count per 
unit flux density at 4.85 GHz, dN/dSs, derived from the PMNS 
catalog, and the distribution of 3 1 GHz to 4.85 GHz flux density 
ratios, 531/55, derived from new observations for this purpose 
with the OVRO 40 m telescope in Ka band (paper in prepa- 
ration). We proceed to calculate the statistical correction for 
unconstrained residual point sources ^31 > 1 mJy using Monte 
Carlo techniques; we generate random point source distribu- 
tions at 4.85 GHz using dN/dSs and statistically extrapolate the 
flux density of each source to each of our ten frequency chan- 
nels using 531/55. These simulated point sources are super- 
imposed on CMB temperature fluctuations and observed with 
DASI simulation software; a power spectrum is then generated 
with the analysis software. The resulting mean amplitudes and 
uncertainties of the residual point source contribution to the 
nine band powers are [20 ±70, 70 ±80, 90 ±70, 180 ±70, 
240 ± 80, 330 ± 100, 400 ± 100, 500 ± 170, 430 ± 170] nK^. 
The reported uncertainties are due to sky sample variance of 
the point source population in the simulations, uncertainty in 
dN/dSs, and uncertainty in 531/55. The residual point source 
contribution diminishes in the ninth band since that band power 
is dominated by visibilities from the highest frequency channels 
where the average point source flux density is lower relative to 
its mean flux density across all ten frequency channels. We use 
these statistically estimated amphtudes and uncertainties to ad- 
just our CMB band-power estimates and uncertainties reported 
below. 

6. RESULTS 

The CMB angular power spectrum from the first season of 
DASI data is shown in Figure 1, with maximum hkehhood es- 
timates of nine band powers, piecewise flat in /(/-I- l)C//(27r), 
spanning the range I =100-900. Adjacent bands are anticor- 
related at the 20% level. In addition, we show an alternate 
analysis of the same data, for nine bands shifted to the right 
with respect to the original band edges, in order to demonstrate 
the robustness of the analysis against possible effects due to the 
anticorrelation of adjacent bands. Note that these two analy- 
ses use the same data to estimate band powers in two differ- 
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ent piecewise-flat theoretical power spectra; only the first nine- 
band analysis (filled circles) is used for the cosmological pa- 
rameter estimation described in Paper III. While increasing the 
number of bands above nine may in principle provide more in- 
formation about the underlying power spectrum, we have found 
that this does not significantly improve our ability to constrain 
cosmological parameters (see Paper HI). 

In a separate analysis, we fit for the maximum likelihood 
value of an additional parameter, the temperature spectral index 
of the fluctuations, 0, where T x . Fitting a single spectral 
index for all nine bands, we find /? = -0.1 ± 0.2 (1 a), while 
fitting a separate spectral index for / < 500 and / > 500 yields 
= -0.2 ± 0.3 and 0.0 ± 0.4 respectively, indicating the fluctu- 
ation power is consistent with CMB. 

Values and marginal uncertainties for the angular power 
spectrum in the primary nine bands are given in Table 1. The 
center and e~'/^ widths of the bands are calculated using band- 
power window functions adapted from Knox (1999) which are 
plotted in Paper III. These are the relevant window functions 
for calculating the expectation value of the band power given 
a theoretical power spectrum. We give the ratio of the un- 
certainty due to sky sample variance to the uncertainty due to 
noise, a^ja^, estimated using the offset log-normal formalism 
of Bond, Jaffe, & Knox (2000). In their notation, as/an is given 
by Cb/xb, where Cb is the band power estimate expressed as 
1{1 + \)Ci/{2-k) and Xb is proportional to the instrument noise 
contribution to the band-power uncertainty. These values may 
be used to estimate the non-Gaussianity in the band-power 
marginal likelihood distributions for parameter estimation cal- 
culations — asymmetric uncertainties due to non-Gaussianity 
are negligible for most of our band powers and we do not plot 
them here. We also tabulate the band-power correlation matrix 
(Table 2). All of the data products necessary for performing 
cosmological parameter estimation from this data are available 
at our website' . 

6.1. Consistency Tests 

We perform three types of tests to check the consistency of 
the data: i) tests on the difference between two visibility 
data vectors constructed from observations of the same fields on 
the sky, ii) construction of a nine-band power spectrum of the 
epoch-differenced visibility data vector, to test for significant 
deviation from zero power, and iii) tests on the difference 
between two power spectra constructed from independent fields 
on the sky. In the second and third types of test, we increase 
the number of frequency channels which are combined in the 
data vector in order to reduce the computational time required 
to produce the power spectra. This increased data compression 
yields a power spectrum similar to the one reported above. 

Of the three types of test, the first is the most powerful 
tool for detecting non-Gaussianity or incorrect estimates of the 
noise. The reduced statistic is 

xVA^ = (Al - ^ifiCni +Cni+Cc,)-' (Al - A2)/Ar, (11) 

where Ai and A2 are the two data vectors, C„i and C„2 are the 
(diagonal) noise covariance matrices, Ccg is the same ground 
constraint matrix that is used in the power spectrum likelihood 
analysis, and = dim A -rank Ccg are the degrees of freedom. 
We spUt the visibilities between the two epochs of available 
observations for each field row, yielding x^/^ = 1 03. This x^ 
value is significant given the A' ~ 3 x 10"* degrees of freedom — 
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Fig. 1 . — The angular power spectrum from the first season of DASI observations, plotted in nine bands (filled circles). We have analyzed the same data in nine 
bands shifted to the right (open circles). The alternate set of band powers are shown to demonstrate the robustness of the likelihood analysis procedure. To extract 
cosmological parameters (Paper III), only the nine bands shown in the primary analysis (filled circles) are used. Adjacent bands are anticorrelated at the 20% level 
(see Table 2). In addition to the uncertainties plotted above, there is a calibration uncertainty of 8%, expressed as a fractional uncertainty on the C; band powers (4% 
in AT/T), which is completely correlated across all bands due to the combined flux scale and beam uncertainties. 



Table 1 

Angular power spectrum band powers and uncertainties 



'eff 




; b 

% 


/(/+l)C//(27r) (pK^) 




/(/+l)Cfi-J^V(27r) (fiKY 


118 


104 


167 


3770 ±820 


23.6 


-250 ±160 


203 


173 


255 


5280 ±550 


31.6 


140 ±120 


289 


261 


342 


3660 ±340 


187 


120 ±120 


377 


342 


418 


1650 ±200 


7.3 


160 ±140 


465 


418 


500 


1890 ±220 


4.3 


70 ±240 


553 


506 


594 


2840 ±290 


4.0 


0±300 


641 


600 


676 


1670 ±270 


2.3 


120 ±420 


725 


676 


757 


2010 ±350 


1.7 


-90 ±580 


837 


763 


864 


2320 ±450 


1.1 


-490 ±850 



The uncertainties listed above do not include flux scale and beam calibration uncertainties. The total combined calibration uncertainty 
is 8% (1 a), expressed as a fractional uncertainty on the C/ band powers (4% in AT/T). 



"/eff is the band-power window function weighted mean multipole moment (see text). 
*/l and /h are the low and high e"'/- points of the band-power window function. 

^(Ts/ct„ is the ratio of the uncertainty attributable to sky sample variance to the uncertainty attributable to noise (see text). 
''Power spectrum of the epoch-differenced data vector described in §6. 1 . 



it indicates that the noise may be slightly non-Gaussian. In fact, 
we see improvement of this statistic if we increase the sever- 
ity of the lunar cuts, but the effect on the power spectrum is 
negligible. It may also indicate that we slightly underestimate 
the noise of the data. However, the uncertainties in all bands 



are dominated by sky sample variance, rather than instrument 
noise, making the power spectrum robust against a noise under- 
estimate of this magnitude. 

A power spectrum in nine bands was created from the epoch- 
differenced data vector, and tested for deviation from zero 
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power using a statistic, with the result x^/N = 9.5/9, which 
is consistent within the 68% confidence interval. The band 
powers for the epoch-differenced power spectrum are given in 
Table 1. 

We use a statistic to test the consistency between power 
spectra generated from each of the four field rows. The x^ 
statistic is constructed from the difference between two power 
spectra which sample independent sky, x^ = C^i - T^2ViPi + 
PiT^iT^] -I?2), where I^i,!?! are the band-power vectors and 
Pi.Pj are the band-power covariance (inverse Fisher) matrices. 
The non-Gaussianity of the DASI power spectrum uncertain- 
ties is small, which justifies using a x^ statistic; we have tested 
its validity with Monte Carlo techniques on simulated data and 
have not found a significant deviation from a x^ distribution. 
The resulting values, with format x^/N (x^ cumulative distri- 
bution function percentile), are: 14.9/9 (91%), 13.3/9 (85%), 
10.9/9 (72%), 4.5/9 (12%), 6.2/9 (28%) and 3.7/9 (7%) for 
the A-B, A-C, A-D, B-C, B-D, C-D differenced field row 
pairs, respectively. The power spectra of the four field rows are 
in reasonable agreement. 

To test the efficacy of the point source constraints described 
in §5, we split the data in each field row between the four fields 
with the highest and four with the lowest detected point source 
flux, and we create two power spectra from the two sets of com- 
bined fields. The x^/^ value for the difference between the two 
power spectra is 11.5/9 (75%) indicating they are consistent 
within the 68% confidence interval. 

Although point sources are the foreground of primary con- 
cern for DASI, constraint matrices are demonstrably effective 
in removing this point source power, and the consistency tests 
above show that the power spectra from sets of fields with very 
different point source flux contributions are in good agreement 
after the constraint matrix is applied. 

6.2. Diffuse Foregrounds 

Diffuse foregrounds are expected to be low in the region of 
the DASI fields at our observing frequency (see Paper 1). To 
check this, we place limits on the contribution of diffuse fore- 
grounds to the power spectrum by creating constraint matrices 
from foreground templates. The constraint matrix formalism 
is a powerful technique for placing limits on foregrounds with 
a known relative intensity distribution, since it allows for arbi- 
trary scaling of the template ampUtude and spectral index, with- 
out knowledge of these quantities at microwave frequencies. 
We create foreground images centered on each of the DASI 



fields from the cleaned IRAS 100 fim maps of Finkbeiner, 
Davis, & Schlegel (1999), cleaned 408 MHz Haslam survey 
maps (Haslam et al. 1981; Finkbeiner 2001), and Ha maps 
(Gaustad et al. 2000; McCuUough 2001). These images are 
converted to visibility template vectors with the DASI simula- 
tion software. We marginahze over modes in the data which 
match the templates using the constraint matrix formalism de- 
scribed in §5. We constrain an arbitrary template amplitude and 
spectral index for each DASI field. With the addition of all of 
these foreground constraints, the maximum change in a band 
power is 3.3%, with most bands changing by less than 1%. 

The Haslam map has a resolution of ~ 1°, making it inad- 
equate as a template for multipole moments > 200, however, 
the power spectrum of synchrotron emission is expected to de- 
crease with / (Tegmark & Efstathiou 1996). Also, the Ha im- 
ages show very low emission in the region of the DASI fields, 
and are of questionable use as a template. As a second method 
of characterizing possible free-free emission, we convert the 
Ha images to brightness temperature at our frequencies assum- 
ing a gas temperature of 10"* K (Kulkarni & Heiles 1988). Sub- 
sequent power spectrum analysis of the resulting image visibil- 
ity templates yields a < 1% contribution to our band powers in 
all bands. 

We conclude that dust, free-free, and synchrotron emission, 
as well as emission with any spectral index that is correlated 
with these templates, such as spinning dust grain emission 
(Draine & Lazarian 1998), make a negligible contribution to 
the CMB power spectrum presented here. 

7. CONCLUSION 

In its initial season, the Degree Angular Scale Interferome- 
ter has successfully measured the angular power spectrum of 
the CMB over the range I = 100-900 in nine bands with high 
precision. The interferometer provides a simple and direct mea- 
surement of the power spectrum, with systematic effects, cali- 
bration methods, and beam uncertainties which are very differ- 
ent from single-dish experiments. We have made extensive use 
of constraint matrices in the analysis as a simple method for 
projecting out undesired signals in the data, including ground- 
signal common modes and point sources with arbitrary spectral 
index. The constraint matrix formalism is also used as a power- 
ful test of correlations with foreground templates with arbitrary 
flux and spectral index scaling; we find no evidence of diffuse 
foregrounds in the data. 

We see strong evidence for both first and second peaks in 



Table 2 

Correlation coefficient matrix for the DASI band powers given in Table 1 



-5.87 X 10-^ 
0.0286 
-0.196 
1 



1.09 X 10"^ -3.01 X 10-* 
-6.65 X 10-3 3.90 X lO'^ 
0.0372 -8.84 x 10"^ 
-0.234 0.0334 
1 -0.193 
1 



-1.08 X 10"* -1.24 X 10-* -1.99x10-* 



-0.243 
1 



0.0349 
-0.182 
1 



-1.00 x 10-3 
-1.61 X 10-3 
-0.0149 
0.0247 
-0.219 
1 



-6.88 X 10-4 
-2.82 X 10-3 
-5.52 X 10-3 
-0.0193 
0.0394 
-0.275 
1 



-4.85 X 10-* 
-2.74 X 10-3 
-6.69 X 10-3 
-8.22 X 10-3 
-0.0200 
0.0339 
-0.220 
1 
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the angular power spectrum at / ~ 200 and I ~ 550, respec- 
tively, and a rise in power at / ~ 800 that is suggestive of a 
third. The detection of harmonic peaks in the power spectrum 
is a resounding confirmation that sub-degree scale anisotropy in 
the CMB is the result of gravitationally driven acoustic oscilla- 
tions such as those which arise naturally in adiabatic inflation- 
ary theories. In addition, the rise in power in the region of the 
predicted third peak strongly supports, from CMB data alone, 
the presence of dark matter in the Universe. 

In this analysis, the authors were guided by the early efforts 
of Martin White, and by numerous enlightening conversations 
with Lloyd Knox, to whom we give our thanks. We are grateful 



for the efforts of Stephan Meyer who, as the director of the 
Center for Astrophysical Research in Antarctica (CARA), lent 
great support to the project. We are indebted to our intrepid 
winterover crew, John Yamasaki and Gene Davidson, for keep- 
ing the telescope in working order, and to Ethan Schartman for 
his extensive hardware contributions. We thank Peter McCul- 
lough and John Gaustad for providing preliminary data from 
the Ha Sky Survey. This research is supported by the National 
Science Foundation under a cooperative agreement (NSF OPP 
89-20223) with CARA, a National Science Foundation Science 
and Technology Center. Support at Caltech is provided by NSF 
grants AST 94-13935 and AST 98-02989. 
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